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ON THE SOLUTION OF A CLASS OF EQUATIONS 
IN QUATERNIONS. 


[Philosophical Magazine, xvii. (1884), pp. 392—-397.] 


THE general equation of the degree w in Quaternions or Binary Matrices 
is obviously w‘, but in certain cases some of these roots evaporate and go off 
to infinity. The only equation considered by Sir William Hamilton in his 
Lectures is the Quadratic Equation of a form which I call unilateral, because 
the quaternion coefficients in it are supposed all to lie on the same side of the 
unknown quantity. I propose here to show how Hamilton’s equation, and 
indeed a unilateral one of any order, may be solved by a general algebraical 
method and the number of its roots determined, 


It will be convenient to begin by setting out certain general equations 
relating to any two binary matrices m, n. 


Writing the determinant of æ + ym+zn under the form 
a + 2bæy + 2cwz + dy? + 2eyz + fe 


(b, c, d, e, f, thus constituting what I call the parameters of the corpus m, n), 
we have universally 


m—2bm+d=0, n?—2cn+f=0, d(m ny - 2e(m—n)+f= 0. 
Moreover if m,n receive the scalar increments p, v; d, e, f become respectively 
d—2ub+ p, e-—po—vb+pr, f—2ve+r* 

Let us begin with Hamilton’s form, say 
a —2pa+q=0, 
and suppose 
| a —2Ba+ D=0, 
where B, D are scalars to be determined, 
_ Let b, c, d, e, f be the five known parameters of the corpus p,q. Then, 
since 
(p — B)7(q— D) = 2a, 
8. IV. 15 
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we shall have [cf. p. 188 above] 


4 (d — 2bB + B?) a? —4(e—bD — cB + BD)a+ f—2D+D?=0. 


Hence, writing B-b=u, D-c=», 
d—lW=a, e—be=B, f—c?=y, 
we have wW+a=rA, wt+R=2rA(ut+b), wP+y=4rA(v+ 0c). 
From the last two equations, eliminating v, there results 
(2Au — 2br— BP — 40 (2Au — 26A— B) u + (y — 4cr) w = 0. 
Hence substituting A — a for u’, 
(4A? + 4cr,— y) (A — a) — (2bà — BP? =0. 
We have thus six values of u, namely 


+ /(A— a) 


(where à has three values), to which correspond six values of v, namely 


4 2r1b— B. 
EATA) 


and, finally, 2¢=(p—u—b)(q—v+0) 
= {(p — b} —w}"(p— b+ u)(q—e—2), 


_ pa—(c+v)p—(b—u)qt(b—u)(ct+»), 
eA oo MBan Iso, eae ee aod ot. nainat 


[32 


which equation gives six values for æ, and shows that ten have evaporated. 


It is easy to account & priori for the solution depending only upon a 


cubic in w?. 


For a?—2p2+q=0 is the same as 7?—2yp+q=0, where y=—a@+ 2p. 
But obviously, from the nature of the process for determining them, B and C 
are independent of the side of the unknown on which the first coefficient lies. 
Hence the actual B will be associated with B’, B’ being what B becomes 


when æ becomes — æ + 2p, which is obviously — B + 2b. 


Hence with any value of B— b, which is u, is associated a corresponding 


B — b, which is — u. 


I will now proceed to apply a similar or the same method to the 
trinomial cubic equation in quaternions (or binary quantity) ° + pæ — q =0, 


with a view to ascertain the number of its roots. 


Retaining the same notation as before, and still supposing 
æ — 2B + D = 0, 
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we obtain æ + (D— 4B?) + 2BD =0, 
_ :¢+2BD , 
and o> ge aE oD 4R?—D 


Hence (4B? — DY} — 2b (4B? — DY) + d} a? 

— 2 {2 (4B? — D) BD —c (4B?— D) — 2bBD + e} æ 

+ 4B2D? —- 2BD + f= 0. | 
Hence we may write 

(4B? — Dy — 2b (4B? — D)+d=2, 
2 (4B? — D) BD — c (4B? — D) —2bBD + e =B, 
| 4BD? — 2cBD +f =)D; 

from which equations B and Dare to be determined. Eliminating between 


the first and second and between the first and third of these equations, we 
obtain two equations, of which the arguments are 


oP aB Ps BD BD BD: Dii 
for the one, 
BI TELD BD D BTB B; 1 
for the other. 
Eliminating D by the Dialytic method between these two equations, we 


shall have (using points to signify unexpressed coefficients) the following 
three linear equations in D?, D, 1, namely: 


-BD? +(- B + &.) D+ (- B’ + &e.) =0, 
- BD? +(- B’ + &c.) D + (. B’ + &e.) = 0, 
- BD? + (- B+ &e.) D +(- B+ &e.) = 0. 


Hence in the final equation B rises to the 15th power; and by com- 
bining any two of the above equations, D is given linearly in terms of B; 
and, finally, œ is known from the equation 


_:(p +D -— 4B? — 2b) (q+ 2BD) 
-= GB Dy—2 (4B — D) +d’ 


æ 
and has 15 values. 


A like process may be extended to a unilateral equation (of the Jerrardian 
form) of any degree, say 2 + ga +r = 0. 


Introducing the auxiliary equation with scalar coefficients as before, namely 
æ— 2B + D=0, 


” may be expressed as a function of q, r, B, D; and the term containing the 


* T use ig and a to signify M—1L and LM— respectively. 
15—2 
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highest power of B in the equation for determining B (of which D is a one- 
valued function), when w = 4, will be found to be the determinant 

By a BRE a B 

+ B® 2 Bi B e 

BR Be BRE 

+p? . B? . Bu . B3* 
and a similar determinant will fix the degree of B in the resolving equation 
for any value of œ. Hence the number of solutions of the unilateral equation 
in quaternions of the Jerrardian form of the degree w is w (2w — 1) or 2%? — ao, 
and the evaporation will accordingly be wt — 20? + w, or 

(@?— w) (+ ow- 1). 

Moreover the same method with a slight addition will serve to de- 
termine the roots of the general unilateral equation in quaternions, the 
number of which will be a cubic function of œ, as I propose to show and to 
give its precise value in some future communication, either in this Journal, 
or at all events in the memoir on Universal Algebra now in the course 


of publication, under the form of lectures, in the American Journal of 
Mathematics t. 


I very much question whether the old method of Hamilton, as taught by 
its most consummate masters, Tait in this country, or the late Prof. Benjamin 
Peirce in America, would be found sufficiently plastic to deal effectually with 
an analytical investigation in quaternions of this degree of complexity, so as 
to lead to the formula for the number of solutions of the unilateral equation 
of the Jerrardian form above given. 


I invite my much esteemed and most capable former colleague and former 
pupil, Dr Story, of the Johns Hopkins, and Prof. Stringham, of the University 
of California, who carry on the traditions of the Harvard School, to put the 
power of the old method as compared with the new to this practical test. 


Postscript.—lIf a — 3pæ + 3qa—r=0, 
(where p, q,r are perfectly general matrices of the second order which satisfy 
the general equations 
g—2bg+d=0, gr+rq—2bqg—2bq+2e=0, 7?-—2hq¢+d,=0, 
pat qp— 2bp —2Bq+ 2e=0, p?’—28p+s8=0, 
pr+rp — 2b,p— 28r + 2e,=0), 


* It may readily be seen that the highest term in the equation for finding B is identical with 
the resultant of 
D*+- 24B?D? + 80B4D? and 4BD*—40B*D? + 64B5D - 64B’, 
that is, will be 2°.3.7.19B*; and that the last term (at all events to the sign près) will be 
b48, which is of 4.34+2.2.4 (that is of 28) dimensions in a, and is therefore codimensional 
(as it ought to be) with B*, 
+ It is given in the Postscript below. 
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and if we write #—2Be+D=0, 
:r +3Dp- BD . 
P?= 3° 3Bp AED 
and I find by perfectly easy and straightforward work that B, D may be 
determined by means of the following equations: 


then 


anog E os wae 2eB + 48B?) = 9), 
AEDE 38D) PTP + (e— aB + 3eD — 63BD) = 3B, 


BD? — 2 (b, + 38D) BD + d, + 6De, + 98D? = Dy. 
The order (by which I mean the number of solutions of this system of 
equations) is readily seen to be the same as that of 
D +-BD+-BD=0 
-BD +- B'D +.B =0; 
that is, is the same as the degree in B of B*(B’)?.R, where R is the 
resultant of 
-D +. B+. Band - D+: BD+- B. 
Hence* the number of solutions is 3 + 10 + 8, that is, is 21. 


Practically, therefore, we have now sufficient data to determine the number 
of solutions of a unilateral equation in quaternions of any order œw; for it is 
morally certain that such number is a rational function of w ; and as it cannot 
but be of a lower order than œt, we have only to determine a cubic function 
of œ whose values for œ =0, 1, 2,3 are 0, 1, 6, 21, which is easily found to 
be œ — w? +w; so that the evaporation is œ — œ’ +œ°—w, that is 

(œ + 1)(@? — o). 

Practically also we can solve (subject to hardly needful verification) the 

number of roots of a unilateral equation of the special form 

L + Got? + Gore t+... + Qo = 0. 

For when 0 = œ, we know the number is w°; and when 0 = 1, the number is 
œ + m?—w; consequently if the second differences of the function of (w, @) 
which expresses the number of roots are constant, the value of this function 
when 0=—1 is œ — wœ? +w, which we have found to be the actual number; 
and consequently, if the second differences are not constant, they must be 
Sometimes positive and sometimes negative, which is in the highest degree 
improbable. Hence in all probability it will be found that the required 
number of solutions in the form supposed is (1 + @) œ? — bo. 


I need hardly add that the nine quantities 2b, 2b,, 28; 2e, 2e,, 2e; d, ò, di, 
which occur in the discussion above given of the general unilateral cubic, or, 
Say, rather the ten quantities obtained by adding on to these wnity, are the 

[* See footnote t p. 197 above. ] 


www.rcin.org.pl 


230 On the Solution of a Class of Equations, ete. [32 


ten coefficients of the determinant to the binary matrix (x + py + qz + rt), 
which of course there is not the slightest difficulty in expressing in terms of 
scalar and vector affections of p, q, r and their combinations, if any one 
chooses to regard them as given in quaternion form. 


Scholium. In what precedes it is very requisite to notice that only 
general cases are considered; and that there are multitudinous others which 
escape the direct applicaticn of this method, and do not conform to the rule 
which assigns the number of solutions. Thus, for example, the equation 
x + px = 0, besides the solutions æ = 0, x =—p, will have two others which will 
require the method of the text to be modified in order to determine. Or take 
the most elementary case of all, the simple equation px =q. If pis not vacuous 
(that is, if its determinant when regarded as a matrix, or its modulus when 
regarded as a quaternion, is finite), there is the one solution e=pq. But 
if p is vacuous, then, unless q is also vacuous, the equation is insoluble. If 
q=0, there will be two solutions; one of them x=0, the other æ = con- 
jugate of p in quaternion Arisinvlogy: or ) 

—d; b a; b 
c= se? when PEG: d 
in the language of matrices. If, p still remaining vacuous, g is vacuous but 
not zero, a further condition must be satisfied, namely, if 


p= © and q= a 
the condition is aò + ad — by —cB=0; 
or if p=a+bi+cj+dk and q=at+fit+yt òk, 
the condition is aa+bB+cy+déd= 0. 


When this condition (besides that of q being vacuous) is satisfied, the equation 
px =q is soluble, and pq becomes finite but indeterminate, containing two 
arbitrary constants *. 


* So in general if p,q be two simply vacuous matrices of any order, the condition that the 
equation pa=q may be soluble, or, in other words, that p~'q (a combination of an ideal with a 
vacuous matrix) may be non-ideal, may be shown to be that the determinant to the matrix 
Ap +uq (where à, u are scalar quantities) shall vanish identically—which (p being supposed 
already to be vacuous) involves just as many additional conditions as there are units in the 
order of the matrix. 
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